QED of excitons with microscopic nonlocality in arbitrary-structured dielectrics 
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We have constructed complete quantum theory for an optical process of excitons with microscopic 
nonlocality (nonlocal susceptibility) originating from their center-of-mass motion. This theory pro- 
vides a practical calculation method for arbitrary-structured nano-to-macro scale dielectrics where 
excitons are weakly confined. We obtain good correspondences with underlying theories, semiclas- 
sical microscopic nonlocal theory, and QED theories for dispersive and absorptive materials with 
local susceptibility. 
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I. INTRODUCTION 

In the conventional theories of optical processes in con- 
densed matters, light has been mainly treated classically 
regardless of whether the matter systems are described 
in quantum mechanical terms (semiclassical theory) or 
classical ones. These theories have successfully explained 
a variety of optical phenomena for the classical light or 
the coherent states of photons. However, there is grow- 
ing interest in the quantum electrodynamics (QED) of 
elementary excitations in condensed matters in order to 
discuss optical processes for nonclassical light such as 
entangled states, single photons, squeezed states, cavity 
photons, and so on. The relevant experiments have al- 
ready been reported, for example, the entangled-photon 
generation via biexcitons (excitonic molecules), 1 trig- 
gered single photon generation from bound excitons in 
a semiconductor^ and the squeezing of cavity polaritons 
in semiconductor microcavities£ The quantization of a 
radiation field has been studied for a long time not only 
in a vacuum^ but also in the medium characterized by a 
frequency- independent dielectric constant. On the other 
hand, Hopfield has systematically discussed the eigen- 
states of exciton-photon systems or exciton-polaritons, 5 
which have a frequency dependence in a susceptibility 
x(w) or in a dielectric function e(u>) — 1 + x( w ) as seen 
from their dispersion relation lo 2 €{lo) — c 2 k 2 . Although 
in his treatment, x(u>) included only a real part, sus- 
ceptibility is generally represented as a complex function 
satisfying the Kramers-Kronig relations. In addition, its 
imaginary part, causing damping effects, cannot be ne- 
glected in the discussion of resonant optical processes of 
elementary excitations in condensed matters. 

The quantization of the electromagnetic fields in such 
dispersive and absorptive dielectrics has been system- 
atically carried out for homogeneous media by Huttncr 
and Barnett (HB)p and for inhomogeneous 3D ones by 
Suttorp and Wubs (SW)^ In the for mer scheme, dis- 
persive dielectrics are described using the Hopfield po- 
lariton model, 5 , i.e., polarizable harmonic oscillators in- 



teracting with a radiation field, and absorption is con- 
sidered using a reservoir of oscillators interacting with 
the polarizable ones. The electromagnetic fields are de- 
scribed in terms of the eigen operators derived from 
the diagonalization of a Hamiltonian. In the expression 
of those fields, there exists a complex dielectric func- 
tion e(u>) represented by system parameters with satis- 
fying the Kramers-Kronig relations. All the quantum 
mechanical properties of the electromagnetic fields are 
characterized by this dielectric function. The pioneer- 
ing work of HB stimulated various theoretical studies 
associated with the QED of dispersive and absorptive 
dielectrics, for example, the spontaneous decayS*? input- 
output relations ; 10 ' 11 ! 12 and quantization in amplifying, 
anisotropic, magnetic, or nonlinear medial On the other 
hand, SW have carried out the quantization of the elec- 
tromagnetic fields in arbitrary-structured 3D dielectrics 
by using the Laplace-transformatin technique, 7,14 which 
is completely different from the quantization scheme of 
HB. Around the same time, the diagonalization of the 
Hamiltonian of SW has been performed by Suttorp and 
van Wonderen^ In these schemes, complex dielectric 
function e(r, tu) depends on spatial position r of a medium 
and radiation frequency u>. 

In the above QED theories and also in semiclassical 
ones, a dielectric function is usually treated as local form 
e(r,u>) with respect to the spatial position. However, 
in general, the optical susceptibility has a nonlocal form 
as x(r, r',ct>), which characterizes polarization P(r,u;) at 
position r induced by electric field E(r',u;) at different 
position r' as 



P(r,a;) = e Q J dr' xO, r', o;)E(r', w). 



(1) 



This microscopic nonlocality originates from the spa- 
tial spreading of the wave function of elementary exci- 
tations or, particularly for excitons in semiconductors, 
their center-of-mass motion with a finite translational 
mass. Usually, such a nonlocality is not considered in the 
discussion of macroscopic materials; this is because the 



2 



coherence length of elementary excitations is generally 
much shorter than the spatial scale of materials. There- 
fore, only the averaged values of physical quantities over 
the coherence volume are reflected in observation, and 
the microscopic nonlocal effect is not apparent. How- 
ever, in order to discuss the excitons in inhomogeneous 
media, we must suppose that the motion of excitons has 
considerably long coherence and microscopic nonlocality 
becomes important even for bulk materials, as explained 
below. 

In the case of homogeneous media, nonlocal suscepti- 
bility depends only on the difference r — r' of the two po- 
sitions; then, Eq. (TTJ) is rewritten in the reciprocal space 
as 

P(k,w) = eox(k,w)E(k,w). (2) 

In this way, susceptibility x(k, ui) has a wavevector de- 
pendence and the w-dependence even for homogeneous 
media when it has a microscopic nonlocality. This k- 
dependence gives more than one propagating or evanes- 
cent modes for a single frequency satisfying 

uj 2 e(k,Lo) = c 2 |k| 2 . (3) 

Now, we consider a single exciton state with finite trans- 
lational mass m ex , transverse exciton energy Et, and 
longitudinal one Et, at respective band edges. Since 
the transverse exciton energy is written as i? cx (k) = 
Et + h 2 \\s.\ 2 /2m ex for wavevector k, we can find two 
propagating polariton modes for frw > Et , and one prop- 
agating mode at the polariton band gap Et < fko < Et 
in addition to an evanescent mode. These multiple po- 
lariton modes do not appear in the Hopfield polariton 
model 5 because the excitons were assumed to have infi- 
nite translational mass. As pointed out by Pekar for the 
first time,— it appears that additional boundary condi- 
tions (ABCs) should be introduced besides the Maxwell 
boundary conditions for the unique connection between 
the polariton modes inside a material and the exter- 
nal ones at the interface between two materials. This 
problem is known as the ABC problem; it arises when 
the translational symmetry of a system is broken due 
to surfaces or interfaces. Since Pekar's work, continued 
studies have revealed that this problem can be resolved 
by considering the microscopic boundary conditions of 
the excitonic center-of-mass motion at interfacesj 17 ' 18 i 19 
Nowadays, in the semiclassical framework, a calculation 
method independent from the notation of ABCs is well 
known as an ABC-free theory 2 ^ or a microscopic non- 
local theory 21 i 22 These theories systematically consider 
the nonlocality of susceptibility \(r^r' and various 
linear and nonlinear phenomena in inhomogeneous ma- 
terials have been discussed using them. In particular, for 
nano-structured materials, where the coherence of the 
center-of-mass motion of excitons is maintained in the 
entire material (weak confinement regime), anomalous 
size dependence of their optical processes has been elu- 
cidated. With regard to nanofihns, the nonlocal theory 



has successfully explained their peculiar spectral struc- 
tures originating from the polariton interference. 19 ' 23 ' 24 
Further, with recent developement of fabrication tech- 
nologies of nano samples, various peculiar effects due to 
long-range coherence are appearing through the interplay 
between the spatial structures of electromagnetic and ex- 
citonic waves, such as the resonant enhancement of a 
nonlinear response^ 5 - the interchnage of quantized states 
due to giant radiative shift, 26 and the ultrafast radiative 
decay with femtosecond order 

From these theoretical and experimental results and 
great interest in the nonclassical states of elementary 
excitations 1 ^! 3 - as mentioned above, it is very attractive 
to discuss them in detail for the sake of applications 
to quantum information technologies. Although some 
studies introduced microscopic nonlocality into the QED 
of dispersive and absorptive dielectrics, there remains a 
problem of how to perform calculations in practical ap- 
plications, as shown in Sec. IIII1 The principal purpose 
of this paper is to construct a QED theory providing a 
practical calculation method for excitons weakly confined 
in arbitrary-structured 3D dielectrics considering radia- 
tive and nonradiative relaxations, which is necessary to 
discuss the effects of, for example, material interfaces, 
excitonic confinement in nano structures, and nonradia- 
tive relaxation processes. In this paper, we merge a mi- 
croscopic nonlocal theor y 21 ' 22 and the quantization tech- 
nique of SWi 7 - 

We explain SW's QED theory for dielectrics with lo- 
cal susceptibility in Sec. [IT] and the previously discussed 
QED theories with microscopic nonlocality in Sec. IIIII 
We show our Hamiltonian in Sec. HVi explain our QED 
theory in Sec. [Vl and discuss the nonradiative relaxation 
of excitons in Sec. lVIl Last, we compare our theory with 
other QED ones with microscopic nonlocality in Sec. IVIII 
and summarize the discussion in Sec. IVIIII We will ex- 
plain only the outline of our theory in these sections 
and redundant calculations are shown in appendices. In 
App-El the second quantization of excitonic polarization 
is discussed. In App.[B] we derive the Hamiltonian from 
the microscopic point of view. In App. [C] we extend 
the Maxwell wave equation discussed by SW so that we 
can consider the excitonic polarization. In App. [D] we 
evaluate the commutation relations of w-Fourier trans- 
formed Heisenberg operators. In App. [El we discuss the 
relationship between the retarded correlation functions 
of excitons and self-standing modes or exciton polariton 
modes. In App.[Fl we verify the equal-time commutation 
relations, which are expected from those of Shrodinger 
operators. 

In this paper, we use MKS units and Coulomb gauge. 



II. QED OF DIELECTRICS WITHOUT 
NONLOCALITY 

In this section, we explain the outline of the discussion 
of SW, 7 the QED theory without microscopic nonlocality. 
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Its Hamiltonian is shown in Eq. (|B9[) of the present paper. 
The quantization scheme is based on the motion equation 
of the electric field E(r, t) whose definition in SW theory 
is 



E(r,t) = -^A(r,t) 



V0 b gM)- 



(4) 



A(r,i) is the vector potential and </>b g (r, t) is the 
Coulomb potential. We write the latter and the dielectric 
function £b g (r, w) with a subscript bg for the description 
in the following sections. Since the Coulomb gauge is 
used in this scheme, the vector potential is a transverse 
field satisfying V ■ A(r) = 0, and the second term of 
Eq. represents the longitudinal field. From the La- 
palcc transform of the Heisenberg equations of the sys- 
tem variables, the Maxwell wave equation for the Fourier 
component of the electric field 



E±(r,^) 



is derived as 



f 

2^ 



dte ±1 " t E(r,t) = {E*(r,Lo)y (5) 



V x V x E+(r,w) 5-e bg (r, w)E + (r, w) = ifi Lo3(r,Lo), 

(6) 

where E + (r,<x>) and E _ (r, lo) are called positive- and 
negative- frequency Fourier components of E(r, t), respec- 
tively. We write the Fourier transformed operator with 
a hat ( * ) in this paper. Wave equation ((6|) has the same 
form as the one that appears in the classical electrody- 
namics except for operator J(r, lo) on the right hand side. 
This operator is called the noise current density and in- 
terpreted as a source of the electromagnetic fields or the 
fluctuation caused by absorption. It plays an essential 
role in the series of QED theories for dispersive and ab- 
sorptive dielectrics. The same kind of operator for homo- 
geneous systems has been derived by HB,? and the one 
for inhomogeneous 3D systems has been phenomenologi- 
cally introduced in Ref. [28|. On the other hand, from the 
Lapalce-transformed motion equations of system vari- 
ables, SW have systematically derived the representa- 
tion of J(r, lo), which is written in terms of the canonical 
variables and momenta of the system at t — (we have 
performed a similar calculation in App. [C]). From the 
commutation relations between them, those for J(r, lo) 
have been derived as 



Here, using Green's tensor G(r, r',Lo) satisfying 



V x V x G(r, r', lo) - ^e bg (r, w)G(r, r', lo) = S(r - r')l, 
we can rewrite Maxwell wave equation © as 



E+(r,w) =wqlo J dr' G(r, r', lo) ■ J(r', lo). 



(9) 



From the commutation relations ([7]) of J(r, lo), those of 
E ± (r, lo) can be derived as 

E+(r, lo), E- (r', to')] = S(lo - co')^-Im[G(r, r', lo)], 

(10a) 

(10b) 



E+(r,u;),E+(r',a/) 



0. 



where we use the equivalence shown in Eq. (1.54) of 
Ref. 0: 

f LO 2 

/ ds -2-Im[e bg (s,w)]G(r,s,w) • G*(s,r',w) 
= Im[G(r,r»] (11) 
and the reciplocity relation: 

G(r,r',u,)={G(r',r,u,)}\ (12) 



Eqs. (|T0|) can be understood by the fact that the Green's 
tensor G(r, r', lo) for the Maxwell wave equation ([6|) iden- 
tifies with the Fourier transform of the retarded correla- 
tion function of the electric field except for the factor 
—Hohw 2 , as discussed in Ref. [291: 

— fi n fkj 2 G(r, r', lo) 

dt gMt-f _ i ')([ E (r, t), E(r', *')]>, (13) 



where 9(t) is the Heaviside step function. In this frame- 
work, all we have to do is find Green's tensor satisfy- 
ing Eq. ([5]) in order to discuss the quantum mechanical 
properties of the electromagnetic fields in dielectrics. All 
the information of the material structure is included in 
dielectric functin eb g (r, lo). The form of Green's tensor 
has already been known for various structures with high 
symmetry^ and also can be numerically calculated for 
arbitrary 3D structures.— 



3{t,lo),{3{v',lo')}^ 
= 5(lo - lo')S(t - r')^^Im[e bg (r, w )]l, (7a) 

TT 



3(r,Lo),3(r',Lo') =0 



(7b) 



where [J, J^] is a 3x3 tensor and its (£, £') element implies 
= x, y, z (the same kind of calculation is 
shown in App. |D]of the present paper). 



III. PREVIOUS QED THEORIES WITH 
MICROSCOPIC NONLOCALITY 

The above theory enables us to discuss the linear opti- 
cal process in arbitrary-structured 3D dielectrics charac- 
terized by dielectric function eb g (r,o;). However, in order 
to discuss the materials with microscopic nonlocality, it 
is necessary to consider more general elementary excita- 
tions that cannot be described by harmonic oscillators of 
the Hopfield model. 
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As a pioneering study on a full-quantum theory with 
microscopic nonlocality, Jenkins and Mukamel have dis- 
cussed molecular crystals in d dimensions (d — 1, 2, 3), 32 
where the relative motion of excitons is localized at a 
single molecule and the center-of-mass moves between 
molecules due to dipole-dipole interaction. While their 
theory concentrates on treating the resonant polariza- 
tion without nonradiative relaxation, recently, micro- 
scopic nonlocality is being introduced into the field quan- 
tization in dispersive and absorptive media , 33 ! 34 ' 35 ' 36 ' 37 ! 38 
and some studies have demonstrated the application 
of their theories for specific structures ! 33 ! 38 Stefano et 
al. discussed excitons with the microscopic nonlocal- 
ity in media with spatial translation symmetry broken 
along one dimension, and they practically calculated the 
spatial and frequency dependences of the vacuum-field 
Actuation in a semiconductor quantum well structure. 33 
Thereafter, they extended their theory to an arbitrary 
3D structure,— and discussed the input-output rela- 
tions in scattering systems^ On the other hand, Bech- 
ler performed the field quantization for the homoge- 
neous systems with nonlocality by using the path-integral 
method^ and Suttorp did for the nonlocal, inhomoge- 
neous, and anisotropic systems by using the diagnaliza- 
tion method.— Most recently, Raabe et al. phenomeno- 
logically discussed the nonlocal systems with both dielec- 
tric and magnetic properties^ They propose the use of 
the dielectric approximation with the surface impedance 
method for the practical application of their theory. 

As seen in the above studies, it is safe to say that 
a consistent framework for the field quantization in di- 
electrics with microscopic nonlocality has already been 
established. Thus, the issue of current importance is 
to establish a general and practical calculation method 
applicable to arbitrary-structured 3D system, which is 
desired for the actual applications of the above frame- 
work though interesting applications have already been 
demonstrated in specific situations by Stefano et al. and 
Raabe et al. The essential task for this purpose is the 
derivation of Green's tensor for the Maxwell wave equa- 
tion with nonlocal susceptibility as seen in Eq. (f59|) of 
the present paper. 

In this paper, we provide a practical calculation 
method for Green's tensor for arbitrary structures by us- 
ing the fact that the nonlocal susceptibility has a separa- 
ble form with respect to two positions as seen in Eq. ({54")) . 
This technique has been developed in the semiclassical 
microscopic nonlocal theory 21 ! 22 We extend this theory 
to be able to consider the quantum mechanical properties 
of electromagnetic fields by using the Laplace transfor- 
mation technique of SWi 7 - In other words, we generalize 
the SW theory to media with microscopic nonlocality for 
providing the practical calculation method. The SW the- 
ory is suitable to extend the nonlocal theory because the 
latter consists of two equatons: the Maxwell wave equa- 
tion, which is just the fundamental one in the SW theory 
as seen above, and the motion equation of excitonic po- 
larization, which can easily be derived by the Laplace 



transformation technique. 

We provide the outline of our theory in the following 
sections, and the detailed explanations, including lengthy 
calculations, are given in Apps.lAltolDlto keep strightfor- 
wardncss of the main part. After that, we will compare 
our theory with some of the above-mentioned previously 
discussed ones. 



IV. HAMILTONIAN 

We describe the dielectric materials with resonant 
contributions from excitons with center-of-mass motion 
and the nonresonant ones with local dielectric func- 
tion Cb g (r, to). This treatment is essential for including 
the consideration of the effects arising from the radia- 
tion mode structures modified by the practical dielectric 
structures (with absorption) such as a substrate, a dielec- 
tric multilayer cavity, photonic crystals, and so on, sur- 
rounding excitonic active structures. We explicitly dis- 
cuss optical and nonradiative-relaxation processes of the 
excitons, and the nonresonant backgrounds are treated 
as the same procedure of SWi 7 " The total Hamiltonian 
discussed in the present paper is written as 



H = H e 



Hi, 



(14) 



where H oia describes the radiation field and background 
dielectric medium, H mat represents excitons and a reser- 
voir of oscillators, and H mt is the interactions between 
H cra and i? m at- H cra is just the complete Hamiltonian 
discussed by SW— , and its representation is shown in 
Eq. (|B9|) of the present paper. 

As a result of the discussion in App. [B] the interaction 
Hamlitonian is represented as 



dr 



Url'Afr) 



1 



iV ex (r)A 2 (r) 



+ J dr <^ bg (r)p ex (r) + - J dr ox (r)/9 cx (r). 

(15) 

I ex (r) is the excitonic current density without radia- 
tion contribution — iV cx (r)A(r), i.e., the complete current 
density is written as J ox (r) = I GX (r) — AT 0X (r)A(r) (see 
App. [A| . Pcx(r) is the excitonic charge density and 



t (r) = J dr' 



Pcx(r') 

47T£o|r — 1 



(16) 



is the Coulomb potential. The first and second terms of 
Eq. (fl~5|) represent the interaction between the radiation 
field and excitons. The third term is the Coulomb in- 
teraction between induced charges of excitons and that 
of the background medium. The last term is the one 
between excitonic charges themselves, and is also consid- 
ered as the dipole-dipole interaction between excitonic 
polarizations, or the exchange interaction between elec- 
trons and hole s 22 ' 39 ! 40 (see App.|B|. Although this term 
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usually belongs to matter Hamiltonian fl ma t , we displace 
it into f/jnt because it can also be considered as the inter- 
action between the longitudinal component of the polar- 
ization and that of the electric field. This treatment will 
give us the motion equation of excitons as a simple form 
as Eq. (|28|) , and will take away our explicit consideration 
of the longitudinal-transverse (LT) splitting of the exci- 
ton eigenenergies, because the last term of Eq. (TT5)) is the 
origin of the LT splitting. 

With regard to excitons, generally, we should describe 
them starting from the basis of electrons and holes inter- 
acting with each other and themselves. However, as long 
as we consider a linear optical process of exciotns under 
weak excitation, it is valid to describe electronic systems 
in terms of excitonic eigenstates. In addition, in order to 
describe a nonradiative relaxation process, we consider 
a reservoir of oscillators interacting with excitons. The 
matter Hamiltonian is written as 

/>oo 

H mat = V fe^&t bfi + V / dn {nn 4(0)^(0) 
+ [K + bl] [g^d^n) + «£(n)d+(n)] }, (17) 

where is the annihilation operator of the excitons in 
eigenstate \l with eigenfrequency u^, which does not in- 
clude the LT splitting because we displace the exchange 
interaction between electrons and holes from fl ma t to 
flint. In this paper, we assume that the center-of-mass 
motion of excitons is confined in finite spaces, and index fi 
represents degrees of freedom of not only the relative mo- 
tion but also the translational one. Instead of evaluating 
the commutation relations of 6 M from its representation 
(|A23p with Fermi's commutation relations of electrons 
and holes, we consider the excitons as pure bosons satis- 
fying 



(18a) 
(18b) 



This approximation is valid under weak excitation. On 
the other hand, in Eq. (jTTJ) , d M (f2) is the annihilation 
operator of the reservoir oscillators with frequency Q in- 
teracting with the excitons in state fi, and 9^(0) is its 
coupling parameter. The oscillators are independent of 
each other and satisfy the commutation relations as 



[d M (n),dj t/ (n')] = W^-flOi 
[d M (n),d M /(n')] = o. 

V. QED OF EXCITONS 



(19a) 
(19b) 



In this section, from motion equations of excitons and 
the Maxwell wave equation for our system derived by 
the quantization technique of SW, 7 we show our QED 
theory for excitons by using the technique of the mi- 
croscopic nonlocal theory developed in the semiclassical 
framework ] 21 ! 22 Since we require a lengthy calculation 



based on the SW's scheme to prove the validity of the 
following discussion, we show only the outline of our the- 
ory and the details are discussed in Apps. IClandlDl We 
do not consider the nonradiative relaxation of excitons in 
this section. It will be discussed in the next section. 

Since we consider excitons and also the background 
medium, the exciton-induced longitudinal field also con- 
tributes to the electric field as 

E(r, t) = ~A(t, t) - V0 bg (r, t) ~ V0 ox (r, f). (20) 

This definition is different from that of SW (Eq. (f4]) of 
the present paper). In addition, instread of Eq. (|6j). 
excitonic polarization density P+.(r, uS) appears in the 
Maxwell wave equation as 



V x V x E + (r,w) - — e bg (r,u;)E + (r,u/) 
= i^ ^Jo(r, w) + Mow 2 P+ (r, uj) 



(21) 



(see App.[C]), where Jo(r, ui) is the same kind of operator 
as J(r, u)) in Eq. ([6]), and it satisfies the same commuta- 
tion relations of Eq. (fTJ) : 



= 8{uj - u')S{r - r')^-Im[e bg (r, w )]l (22a) 



J (r,w), J (r',u/) 



0. 



(22b) 



(see App. [P|l . This operator is interpreted as the source 
of the background electromagnetic fields and its defi- 
nition is shown in Eq. (|C32|) . Using Green's tensor 
G(r, r',oj) satisfying Eq. ©, we can rewrite Eq. (|2"Tj) as 

E+(r, u) = E+(r, u) + ^ 2 J dr' G(r, r', w) ■ P+ (r', w), 

(23) 



where 

E+(r, W ) 



dr'G(r,r',w).Jo(r',w) (24) 



is the background electric field satisfying the commuta- 
tion relations 



fe+(r,w),Eo(r>') 



E+(r, W ),E+(r',c'; 



0. 



(25a) 
(25b) 



These are equivalent to relations (Ti"Tj|) of the electric field 
in a local system. This equivalence indicates that the 
behavior of the background field in our system is exactly 
the same as that of the complete electric field in the local 
system. This natural result can be systematically derived 
from the representation of Jo(r, u>) (|C32I) and commuta- 
tion relations of the system variables by using the Laplace 
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transfromation technique. We have verified this equiva- 
lence for the linear optical process of excitons even with 
the nonradiative relaxation of excitons (see App. |D|) . 

Next, we discuss the motion of unknown variable 
P+.(r,(x>) in the Maxwell wave equation (|2ip . The 
second-quantized polarization density is written in terms 
of exciton operator set {6^} as 



Pex(r) 



b U 



where the expansion coefficient ^(r) is 



TJt) 



(26) 



(27) 



(see App. IA"]) . is the transition dipole moment, e M 
is a unit vector in the polarization direction, and G M (r) 
is the wave function of the center-of-mass motion in ex- 
citon state /x. Since we assume the weak confinement 
regime, Vn approximately depends only on the rela- 
tive motion of excitons and is related with LT splitting 
^lt = l^lV^bg- Neglecting the reservoir oscillators 
and assumimg w ~ w^,, the Fourier-transformed Heisen- 
berg equation of excitons is derived from Hamiltonian 
(fl5l) and Q7H as 



{Hu„ -tiu- IS) S M (w) = J dr 7»; (r) • E+ (r, u) (28) 

(see Add. ID lj) . where infinitesimal damping id is added to 
match the discussion with the nonradiative relaxation in 
the next section. Under the rotating wave approximation 
(RWA) , the positive- frequency Fourier component of the 
polarization density (126|) is written as 



P+ x (r, W )=^P /i (r)S /J (o 



(29) 



Substituting Eq. (|28p into (|29| . we obtain the nonlocal 
form of the polarization density as 

P+(r,a,) = e / dr'x cx (r,r»-E+(r',^), (30) 



(31) 



where the susceptibility tensor is defined as 



Xex(r,r',w) = — V 



hwn — ftw — i6 



The spatial spreading of the exciton state, the origin of 
the nonlocality, is reflected through polarization coeffi- 
cient PjuO") or center-of-mass wave function G M (r). On 
the other hand, the spatial structure of the background 
dielectrics is characterized by dielectric function eb g (r, ui) 
in Maxwell wave equation (|21[) and in commutation re- 
lations (|2"2")l . In our framework, we can discuss arbitrary- 
structured exciton motions and background dielectrics as 
seen below. 

In order to discuss the optical process of excitons, we 
must simultaneously solve Maxwell wave equation (12ip 



and the motion equation of polarization density (130)) 
to determine unknown variables E + (r,w) and Pj^r, ui). 
Substituting Eq. (f3"0")) into (|21[) . we obtain the nonlocal 
wave equation as 

uo 2 

V x V x E+(r,w) - — e bg (r,w)E+(r,c i ;) 
c z 

dr' Xox(r, r', ui) ■ E+(r', ui) = ut ^Jo(r, ui). 

(32) 

This is the same equation discussed by Savasta et a h 36 i 37 
and also has the same form as that of Raabe et al.— 
However, it appears very difficult to solve this nonlocal 
equation, and there remains a problem to derive Green's 
tensor for this equation in the practical application of 
their theories. This problem can be solved by using the 
fact that nonlocal susceptibility |3T|) has a separable form 
with respect to r and r'. One scheme is to directly de- 
rive Green's tensor for Eq. (|32|) as discussed in Ref. [HI . 
and the other is to reduce this nonlocal problem into a 
simultaneous linear equation set i 21 ' 22 In our QED theory, 
we adopt the latter scheme becuase it provides not only 
Green's tensor of the former but also much interesting 
information on exciton-polariton systems. 

Substituting the representation of the electric field (j23|) 
into the motion equatin of excitons (|28p by expanding 
P+^i-jw) as Eq. (f2"9")) , we obtain the linear equation set 
determing exciton amplitudes {b^(uj)} as 

y"l [(H -fa- i<5)<W' + ^W' M] V M 



{Ttw^-fvu)-iS)Vf>{u)) 



(33) 



This has the same form of the self-consistent equation set 
in the semiclassical nonlocal theoryi 21 i 22 The operator on 
the right hand side 



1 



fan — fa — lS 



dr7>*(r).E+(r,oO (34) 



has the same form of Eq. (f2"B"|) replacing E + (r,ijj) with 
Eq (r,cj); then, it can be interpreted as the exciton am- 
plitude directly induced by the background electric field. 
Here, we use the word "directly" to mean that b l °\uj) 
does not include the diffusion of the exciton amplitudes 
via the electromagnetic fields. Such effect is reflected in 
the correction term of eigenenergy hw^ in Eq. (|33|) : 



J dr y"dr' p;(r)-G(r,r»-7V(r') 



(35) 

This value represents the exciton-exciton interaction via 
the electromagnetic fields, i.e., the polarization at r' in- 
duces electric field, and later it induces another polariza- 
tion at r. The interaction between the transverse fields 
is the retarded interaction, and the one between the lon- 
gitudinal fields is interpreted as the Coulomb interaction 
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between induced charges. The latter is just the exchange 
interaction between electrons and holes, which we dis- 
place from -Hmat to H- mt , and it gives the LT splitting 
of the exciton eigenenergies. The radiative relaxation of 
excitons and the transition between exciton states are 
described by the imaginary part of (uj). 

Writing the coefficient matrix element of Eq. (|33|) as 



the self-consistent equation set is rewritten as 

]TS^MVH = / dr 7>*(r) -E+(r, W ). (37) 

By using the inverse matrix W(u;) = [S(oj)] _1 with the 
basis of exciton states, we obtain the representation of 
exciton operators as 

&» = X WW M / dr ^ (r) • E+ (r, w). (38) 

We can describe all the physical variables in terms of 
these operators and E^(r, u>). For example, the excitonic 



polarization is written as Eq. (|29ll and the electric field 
(|23| as 

E+(r, u) = E+(r, u) + J2 wji^w), (39) 

A" 

where the coefficients are defined as 

£ M (r,u,)=/W / dr'G(r,r', W ).P M (r'), (40) 

F^v,uj) = ^ 2 J dr'G(r,r',a;)-p;(r'). (41) 

The latter will appear in Eq. The exciton oper- 

ators are represented by the background electric field 
Eg (r, w), whose commutation relations are described by 
Green's tensor G(r,r',w) as shown in Eq. (|25|) . Since 
there is no problem in deriving Green's tensor as ex- 
plained above, the additional works are to perform inte- 
grations (|35|) , (|40|) , and (|4T|) , and to derive inverse matrix 
W(uj) = [S(w)] _1 . These calculations can be performed 
straightforwardly. Therefore, in our framework, for a 
given wave function of excitons' center-of-mass motion 
{G AI (r)} and background dielectric function eb g (r, uj), we 
can discuss the QED of excitonic materials with the mi- 
croscopic nonlocality. 

When we do not consider the nonradiative relaxation 
of excitons, all the quantum mechanical properties are 
described by commutation relations (|22[) of noise cur- 
rent density Jo(r, w). Using relations (|25|) derived from 
Eqs. (|2"2"|) . the commutation relations of exciton operators 
(1551 are evaluated as 



5^-u>') — [w^^)-w;,^)], 



6 M (w),V(w') =0 



(42a) 
(42b) 



This result indicates that the elements of inverse matrix 
\N(ui) of the self-consistent equation set identify with the 
Fourier transforms of the retarded correlation functions 
of excitons except for the factor —H (see App.|E|. In addi- 
tion, we can obtain the commutation relations of electric 
field operators (|39jl : 



E + (r,w),E-(r',a/) = 5(u - u') 
[i^fa 2 



i2ir 

E+(r jW ),E+(rV) 



[G re „(r, r', u) - {G rcn (r', r, w)}**] , (43a) 
0, (43b) 



where G rcn (r, r', uj) is defined as 
G ren (r,r',w) = G(r, r',u) 

±-J2 £»{*,") ^WM *V(r»- (44) 



We can find that if this tensor satisfies the equation 

uj 2 

V x V x G ren (r, r', w) -e bg (r, w)G ron (r, r', w) 

c z 

" ^ / dr" Xc X (r, r", w ) • G rc „(r", r', w ) = <5(r - r')l, 

(45) 

then it can be interpreted as Green's tensor for nonlocal 
Maxwell wave equation (|32|) . This tensor, which renor- 
malizes the linear optical process of excitons with the mi- 
croscopic nonlocality, is also calculated directly from the 
above nonlocal equation.— The breakdown of reciprocity 
relation (|T2|) for G ron (r, r', uj) arises from anisotropic sus- 
ceptibility tensor (|3 If) of the excitonic polarization. By 
using G ren (r, r', uj), electric field operator ([39]) is also 
written as 

E+(r, uj) = ifx uj / dr' G rcn (r, r', uj) • J (r', uj). (46) 



Comparing Eqs. ©, (HUJ), (g3J), and gBJ, we can find a 
good correspondence between our theory and the QED 
theories for local systems. In addition, we also verify that 
commutation relations (|42p and (|4"31 satisfy the equal- 
time commutation relations expected in the Shrodinger 
representation (see App. |F|) . 



VI. WITH NONRADIATIVE RELAXATION 

Next, we discuss the modulation of the above thoery by 
considering the nonradiative relaxation of excitons. The 
detailed caluclation is shown in App. ID 3[ and the result 
is as follows. The nonradiative relaxation slightly mod- 
ulates the motion equation of excitons from the original 
one (gll) : 

[H - fa - i7 M (w)/2] bf,(uj) 

dr?;(r).E+(r, U ) + ^( W ). (47) 
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7/i (^0 is the nonradiative relaxation width defined as 
Eq. (|D48p in terms of exciton-reservoir interaction co- 
efficient g^iQ,). Operator 2? M (a;) represents the fluctua- 
tion caused by the reservoir oscillators. It is defined as 
Eq. (|D51|) and satisfies the commutation relations 



w'Siu- u')— 7 M (w), (48a) 

Z7T 



(48b) 



This is another source operator of our system and is in- 
dependent from noise current density J (r, w) as 

J (r, = fe(w),{Jo(r, - 0. (49) 



Substituting representation (|23f of the electric field into 
Eq. Wf\ . we obtain the self-consistent equation set as 



dr K(r).E+(r,c)+^H, 



(50) 



'M- (51) 



where the coefficient matrix element is 

= - fiw- i7/*M/2] + A 

The above result indicates that we can easily intro- 
duce the nonradiative relaxation into the discussion of 
the previous section. All we have to do is introduce re- 
laxation width Jfi(u>) into coefficient matrix (|5ip of the 
self-consistent equation set. By deriving its inverse ma- 
trix w abs (w) = [s abs H]-\ we obtain the commutation 
relations for excitons as 



S»,{VK)} f 



S(u </) — [Wffi M - {W^(u)}*] , (52a) 



fyuM, v (<*/) 



0. 



(52b) 



These have the same forms as those of Eqs. (|42p . 

On the other hand, substituting Eq. (|47j) into Eq. ([29]). 
the excitonic polarization (|30p is rewritten as 



P+ (r,w) = £o / dr' X c abs (r,r',c) • E+(r,w) 

nW^/i(w) 



Tujp — huj — i7 M (u>)/2' 



(53) 



where the nonlocal susceptibility is also rewritten as 

,abs/„ „/, A _iy W) 



Xe a x S (r,r', W ) 



(54) 



Substituting Eq. ([55]) into Maxwell wave equation (|2"Tj) . 
nonlocal wave equation (|3"2")) becomes 



V x V x E+(r,w) - ^-e bg (r,w)E + (r,t 



dr' xS 8 (r,r» ■ E+(r» = i^J abs (r, W ), 



(55) 



where we define a new noise current density as 

j abs (r, U ) ^ Jo(r, «) - L, V y ^^ Wo . 
V ; V ; ^ /^-/^-i 7/ ,(w)/2 

(56) 

From commutation relations (|2"2")l , (|4"5| , and (|4"9")l , we can 

obtain the commutation relations for jQ bs (r,(jj) as 

j abs (r,^),{J ab VV)} f 



8{w-u/)- 



j abs (r,^),J abs (r>') 



7T 

abs/„/ , A 



Im[e abs (r, r',uj) 



0, 



(57a) 
(57b) 



where e abs (r, r', w) = l<5(r — r')eb g (r, to) + Xcx S ( r i Y ' ■> w ) 
is the nonlocal dielectric tensor. This is a natural re- 
sult from the fluctuation theorem as discussed in Ref . 1361 
In addition, we also obtain the commutation relatons 
for the electric field as having the same form as that of 
Eq. (|4"3")) by replacing G rcn (r, r', uj) with 



G ab :(r,r',c) = G(r,r» 



Moo; 



(58) 



which satisfies the nonlocal wave equation 

V x V x G ab n s (r,r» - ^e bg (r, c,)G ab n s (r, r', u>) 

- ^ | dr" X c abs (r,r'» • G ab n s (r", r', w ) = <5(r - r')l 

(59) 

and Green's tensor required for the practical application 
of nonlocal QED theories i 3 -^ 



VII. DISCUSSION 

In this paper, we have constructed a QED theory 
for excitons weakly confined in arbitrary-structured di- 
electrics with microscopic nonlocality and nonradiative 
relaxation of excitons. On the other hand, as mentioned 
in Sec. lIIIi the QED of media with the microscopic nonlo- 
cality has already been discussed in a few studies. From 
the viewpoint of a practical application, we compare our 
theory with the studies of Stefano et a h 33 ' 36 i 37 and Raabe 
et slM 

Stefano et al. have discussed quantum-well structures 
of the dispersive and absorptive dielectrics with micro- 
scopic nonlocality in Ref. l33l . and their theory is gener- 
alized to enable the consideration of arbitrary structures 
in Ref. [H and[37]. However, there still remains a problem 
in deriving Green's tensors for the nonlocal wave equa- 
tion as shown in Eq. (|45[) or (|59p in the present paper. 
Our thoery provides a solution to this problem by giving 
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a definite calculation method using Green's tensor sat- 
isfying wave equation (jSJ) for local media and the fact 
that the nonlocal susceptibility has a separable form as 
shown in Eq. {31]) or ((SJ). The problem of Ref. H3 can 
be solved by using our theory because we derive Green's 
tensor (|44|) or ([58]) for arbitrary-structured excitonic po- 
larization and background dielectrics. 

On the other hand, Raabe et al. propose the use of 
the dielectric approximation with the surface impedance 
method for the practical calculation of Green's tensor for 
the nonlocal Maxwell wave equation. In the dielectric 
apprximation, the characteristic length of spatial disper- 
sion (the spatial spreading of the center-of-mass motion 
of excitons) is assumed to be small compared to the spa- 
tial length of materials, and the information outside a 
focusing region is compressed to integrations of the elec- 
tromagnetic fields at the interfaces. Green's tensor can 
be derived using the surface impedance method for a 
given surface impedance or admittance, which just in- 
cludes the outside information. In contrast, our thoery 
provides Green's tensor for given £b g (r, u;) and center- 
of-mass wave functions of excitons by applying only the 
RWA. Since nonlocality becomes essentail only under the 
resonance conditions, the RWA does not impose any sig- 
nificant restriction on our thoery for the discussion of 
nonlocal systems. In addition, our calculaton method 
would be generalized to the one without the RWA as 
performed in the semiclassical nonlocal theor y 21 i 22 

As mentioned in Sec. HI there is grwoing interest in 
the QED of elementary excitations in condensed matters. 
For example, theoretical studies on entangled-photon 
generation via biexcitons have already been performed 
by Savasta et al . 42 i 43 (though the microscopic nonlocal- 
ity was not sufficiently considered in these calculations) 
and also by us by extending our QED theory reported 
in the present paper4i When we discuss such nonlin- 
ear processes of excitons, we must self-consistently treat 
their nonlinear motion equation and the Maxwell wave 
equation. Based on self-consistent equation set ([57)1 or 
(|50| as discussed in the present paper, the new task is to 
solve the equation set with nonlinear terms originating 
from nonlinear processes. On the other hand, based on 
Maxwell wave equation ([52")) or ([55]) with the nonlocal sus- 
ceptibility as discussed in the previously discussed QED 
theorieS ) 33 i 34 ' 35 i 36 ' 37 i 38 we must solve the wave equation 
with both nonlinear and nonlocal susceptibility. Both ap- 
proaches can be performed by applying some techniques 
such as successive approximation; however, they gener- 
ally require much hard work. In such a case, more de- 
tailed and systematic calculation should be performed by 
using the Feynman diagram technique with the correla- 
tion functions of excitons derived using our QED theory. 
In this sense, our scheme will be a powerfull tool to dis- 
cuss the nonlinear quantum optics in condensed matters 
with the microscopic nonlocality. 



VIII. SUMMARY 

By using the quantization technique of Suttorp et 
al.fiii 4 - we have extended the microscopic nonlocal 
theor y 21 i 22 developed in the semiclassical framework to 
enable the consideration of quantum mechanical prop- 
erties of the electromagnetic fields. This theory keeps 
good correspondences with both the nonlocal theory and 
the QED theories for dispersive and absorptive materials 
with local susceptibility. While microscopic nonlocality 
has been considered in the form of nonlocal susceptibility 
in some of the QED theories ^i 3 ^^ 3 - 7 -^ we reduce the 
problem to a linear equation set as discussed in the semi- 
classical frameworki 21 ' 22 In addition, this theory provides 
Green's tensor for the Maxwell wave equation with a non- 
local susceptibility, which is just required for a practical 
application of previous studies. By using our theory, we 
can discuss the QED of excitons weakly confined in nano- 
structures, which are known to show anomalous nonlin- 
ear optical phenomenaj 25 i 45 ' 46 i 47 ' 48 Although this paper 
is devoted to only the linear optical process of excitons, 
we can phenomenologically extend our theory to describe 
nonlinear processes, such as the entangled-photon gener- 
ation via biexcitons in nano-structures^ 4 - In addition, our 
theory has a potential to systematically discuss higher or- 
der nonlinear processes of elementary excitations in con- 
densed matters by using the Feynman diagram technique 
with correlation functions of excitons derived in our cal- 
culation. Based on our QED theory, we are going to 
discuss various optical phenomena which cannot be dis- 
cussed in the semiclassical framework. 
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APPENDIX A: THE SECOND QUANTIZATION 
OF EXCITONIC POLARIZATION 

In this appendix, we provide microscopic descriptions 
of the current density, charge density, and polarization 
density of charged particles. Then, we expand them in 
terms of the electron or exciton operators. We write the 
second-quantized operators with a hat ( * ) in this ap- 
pendix. 

Considering the charged particles with mass mj and 
charge qi at position rj, current density J cp (r) and charge 
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density p cp (r) are written as 

J cp (r) ee£| [US(t - Ti ) + S(r - ri )ri] , (Al) 



Pcp(r) =^gi<y(r-ri) 



(A2) 



Here, due to the interaction with the radiation field (see 
App. [Bj), the momentum of the charged particles is writ- 
ten as 

Pi = m l v l + g»A(r<). (A3) 

Then, current density (|A1 1) includes a contribution from 
the radiation field. In order to expand it in terms of 
electron or exciton operators, we extract the radiation 
contribution from J cp (r): 

I cp (r) = V [ Pi 6(r - ri ) + S(r - r,)Pi] . (A4) 
* — ' 2m,i 

i 

By writing the coefficient of the radiation contribution 

as 

mi 



(A5) 



we can write complete current density (|A1[) as 

J cp (r) =I cp (r)-JV cp (r)A(r). (A6) 

This subtraction of the radiation contribution is dis- 
cussed in section 2.2 of Ref. [22j. 

Next, we expand the above variables in terms of the 
electron operator a v and its wave function Lp 71 (v). The 
field operator is written as 

%j}(r) = ^2a c ip c (r) +^a„</J„(r), (A7) 

C V 

where labels c and v represent the degrees of freedom of 
conduction and valence electrons, respectively. Assuming 
optical excitation of electron-hole pairs, we obtain the 
second-quantized form of the above variables as 

Iex(r) = -r^ V ala c [p*(r)p<p c (r) - <p c (r)p9?*(r)] 

+ H.c, (A8) 

A> cx (r) = t*L V ata c ^(r)^ c (r) + H.c, (A9) 
m ' 

/5 0X (r) = (-e) ^ aja c ^(r)^ c (r) + H.c. (A10) 



These operators are also written in terms of exciton op- 
erators {b^} as polarization density (|26|) : 

icx(r) =5^T M (r)6 M + H.c., (All) 
iV ox (r) = ^ JV„(r)S„ + H.c, (A12) 
p cx (r) = ^ p M (r)S„ + H.c. (A13) 



Instead of evaluating the expansion coefficient of each 
operator, we describe them in terms of 7- , /J (r), the coef- 
ficient of polarization density (|26| . From the relations 



Pex(r) 



-V-Pe X (r), 



Jcx(r) = ^Pcx(r) = ^ Pex(r),i? 
ot in 



(A14) 
(A15) 



and considering weak exciton-photon interaction, i.e., 
H ~ -ffmat and J ex (r) ~ lex(r), we can write the above 
coefficients as 



2T M (r) = -ko^V^r), 
P/iW = -V-P^(r), 



(A16) 
(A17) 
(A18) 



where cj m is the eigenfrequency of excitons. Using above 
operators (|A1 ![) - i| A13[) . the excitonic current density and 
Coulomb potential are respectively written as 



Jex(r) = Iex(r) - iV ex (r)A(r), 
Pcx(r') 



&x(r) = / dr' 



47re |r 



(A19) 
(A20) 



In order to evaluate coefficients (|A16|) - (|A18|) . we derive 
the representation of "P ^ (r) from the microscopic descrip- 
tion of the polarization density. Averaging the polariza- 
tion at lattice point Ro over a unit cell, the microscopic 
description is written as 

P cp (R ) = i I dr^g.r ^(Ro + r-r,), (A21) 



where the integration is over the unit cell and f2 is its vol- 
ume. Explicitly indicating the lattice point of the elec- 
tron states as (?7, R) and assuming their wave function as 
the wannier function w n (r — R), we obtain the second- 
quantized form of the polarization density as 

P OX (R ) = % iRo+R ,a Cl R 0+ R +R / 

c,u,R,R' 



x — ( drw*(r- R') (— e)r w c (r 



R-R) 



We expand this in terms of the exciton operators 



'/>-"< £^ '''».<. >•.!(. ''k.R, ",-.H. H H- 

c,v,Ro,R 



G* 



I- H.c. 

(A22) 

(A23) 



where /i and "J^c^.r are, respectively, the quantum num- 
ber and the wave function of the relative motion of exci- 
tons, and m and G mj R are those of the center-of-mass 
motion. From the completeness of the wave functions, 
we can rewrite the electron-hole operator set as 

%,R ^c : R +R. — ®fi,c,v,n. G m> R b^m. (A24) 
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Using this relation, we can expand Eq. (|A22[) in terms of 
the exciton operators as 



P CX (R ) = ?Vm(R.o) V™ + H - C - 



(A25) 



where the expansion coefficient is written as 

m,Ro+R' ^/a,c,u 5 R 

c,u,R,R' 

x 1 f dr w* v (r - R) (-e)r w c (r - R- R ). (A26) 

Supposing that the spatial variation of the center-of- 
mass wave function is negligible within the extent of 
the electron-hole relative wave function, we can consider 
G m ,R 0+ R' ~ G m ,R . Further, by expanding the integra- 
tion range to the entire crystal region with iterating R', 
we obtain 

c,t>,R, 

x i Jdrw* v (r) (-e)r io c (r - R). (A27) 

Here, if we assume the wave functions to be smooth 
with respect to the spatial position, i.e., G m (Ro) = 

and $ M)C ,u(R) = $^,c,«,r/ V^, w e obta in the 
expansion coefficient of polarization density (|A25[) as 



'P/j.miRo) — 'P/i G m (Ro), 



(A28) 



where 



V» = $ M,c,v(R) / dr u£(r) (-e)r u; c (r - R) 



(A29) 

is the transition dipole moment of exciton band p, and 
its absolute value is related with the LT splitting of the 
exciton eigenenergy as A^ T = \V / e eb g - 



APPENDIX B: DERIVATION OF HAMILTONIAN 

As a model of the background system, i.e., radiation 
field and local dielectrics, we adopt the system discussed 
by SW. 7 Considering the charged particles of App.[2] the 
total Lagrangian is written as 



position: 



- (0bg + 0cp)(Pbg 

l f°° 

+ -p I duo Y 2 



duj u,,X ■ Y„ 



\p~X- 2 - t^ 2 x 2 



Pep J 



-aX + J r 



duo oj 2 Y u 



(B2) 



We omit the descriptions of position dependences. E = 
—A V^bg ~ V0 cp is the electric field, and B = V x A is 
the magnetic induction. X(r) is the amplitude of polar- 
izable harmonic oscillators with density p(r) and eigen- 
frequency u>o(r). These oscillators describe the back- 
ground medium in our QED theory. The polarization 
density, charge density, and current density of the back- 
ground are respectively written as -aX, pbg = V ■ (aX), 
and -aX with position-dependent coefficient a(r). The 
background Coulomb potential is written as 



V = / dr' 



47re r 



47re r - r' 



0b g is related with the longitudinal component of the 
polarization as 



V^bg = 



1 



JL • 



and it also satisfies the Poisson equation 
V72-. __Pbg_ 1 



£0 



-V ■ (aX). 



(B4) 



(B5) 



The damping in the background system is described by 
a reservoir of oscillators interacting with the polarizable 
ones. Y w (r) is the amplitude of the oscillators with fre- 
quency uj, and v w (t) represents the coupling strength. 

From Lagrangian (|B1[) . the canonical momenta of the 
above variables are derived as 



n 



Qu 



Since A and V 



dL . 
— =- = eoA, 
OA 

dL_ 

ax 

dL ■ 
— i — = pY, 
dY u 

dL 

— = rrnvi 
dvi 



pX — aA, 



3»A(rj). 



(B6a) 
(B6b) 
(B6c) 

(B6d) 



cp ) are perpenticular to each 



L = E 



im.r 2 - V{vi) 



other, from the Poisson equation (|B5[) and V 2 
dr C, (Bl) — p cp /eo, the first term of Eq. (|B2[) is rewritten as 



'€]> 



where V(ri) is the one-body potential of the particles 
and C is the Lagrangian density depending on the spatial 



dr ^E 2 

2 



dr 



n 2 i 

+ 2^ b s + ^c P )(pbg + Pbg) 



(B7) 
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After a straightforward calculation, we obtain the Hamil- 
tonian as 



H = H e 



E 



^- { Pi - %A( ri )} 2 + V(r z ) 



dr 



1 



-0cpPcp + 0bgPc 



(B8) 



where 7? cm describes the complete Hamiltonian discussed 
by SW, representing the radiation field and background 
dielectrics with local susceptibility: 



H P . 



dr 



tt2 1 p2 

iL + _L( V x a) 2 + 

2e 2^ V ' 2p 



+/„ -^ Y - 2 

^A 2 + fd W ^X-Q w + 1 
2p Jo P 



2 " 

a. 



-P A 



(B9) 



The first two terms represent the radiation energy, the 
third term is the kinetic energy of the oscillators, and the 
fourth is the potential. The seventh and eighth terms rep- 
resent the interaction between the oscillators and the ra- 
diation field. The eigenfrequency of the oscillators shown 
in the forth term of Eq. (|B9[) is modified as 



LJQ 2 = LOQ 2 + 



PJo 



(BIO) 



by the interaction with the reservoir oscillators, which is 
described as the ninth term. The energy of reservoir is 
the fifth and sixth terms. The last term is the Coulomb 
interaction between the induced charges of backgrounds. 

The kinetic energy of the charged particles, the second 
term of Eq. (|B8j) , is expanded with the expression (|B6dj) 
of their momentum as 



\ " 



1 



— E o m , p 



[Pi - q l A(r i )Y 



4^ 2m, 



[pi ■ A(r t ) + A(r f ) ■ p 4 ] 



from the exciton-associated components of the last three 
terms of Eq. (|B12[) . i.e., the terms proportional to a\a c 
or a\a v but not to aj,a c / or ala v > , which are negligible 
under the weak excitation regime. On the other hand, 
as mentioned in Sec. lIVi we put the exchange interaction 
between electrons and holes into H ln t. It is obtained by 
expanding the fourth term of (1B12|I : 



~ J dr C pp cp 



^ a\a v a\,a c ' / drdr' 



c.v,c ,tr 



, eVc( r )^( r )^'( r ')Vo'(r') 



47reo|r — r'| 



(B13) 



When we use the RWA and assume commutation rela- 
tions (fT5|) of the exciton operators, we can find that the 
Coulomb interaction between the excitonic charges them- 
selves 



- /drdr' P«*( T )P°*P) 



47reo r — r' 



(B14) 



gives exchange expression (IB13I) with a constant energy 
term by expanding p cx (r) using Eq. (|A10[) . While all 
the other terms (• • • ) in Eq. (|B13|) and the second and 
third terms of Eq. (|B12|) belong to H ma t, instead of dis- 
cussing them in detail, we treat the matter Hamiltonian 
as Eq. (fTTj) for a simple description of the linear optical 
process of excitons with nonradiative relaxation. 



APPENDIX C: EXTENSION OF MAXWELL 
WAVE EQUATION 

Here, we extend the Maxwell wave equation discussed 
by SW to enable the consideration of the exciton-induced 
polarization with nonlocal susceptibility. We derive the 
Heisenberg equations of the system variables in Sec. IC 1[ 
calculate their Laplace-transform in Sec. IC 21 and pro- 
vide the Maxwell wave equation for the electric field in 
Sec. EH 



^ 2m 



■A 2 ( ri ). 



(Bll) 



1. Heisenberg equations 



The first term is the kinetic energy without the radiation 
contribution, and the other terms represent the interac- 
tion between the charged particles and the radiation field. 
Here, using the variables defined in Eqs. (|A4[) and (IA5|) . 
we can rewrite Hamiltonian (IB8I) as 



H = H C 



E 



^-p l 2 + V{v l ) 
2m, 



-'cpPcp 



J dr b gPcp - J dr 



Icp ■ A. 2"^ G P"^" 



(B12) 



Expanding these terms with field operator (|A7|1 . we ob- 
tain the first tree terms of interaction Hamiltonian (|15|) 



We derive the Heisenberg equations of the system vari- 
ables and momenta from background Hamiltonian (|B9|) 
and interaction terms (|15[) . The commutation relations 
of the variables are 



[A(r),II(r')] 
[X(r),P(r')] 
[Y w (r),Q w ,(r')] 



where 



6 T (r 



ih 5 T (r-r'), 
ih S(r - r') 1, 
ih 5(u-w') tf(r-r') 1, 



1 S(r 



V'V 



4?r r - 



(CI) 
(C2) 
(C3) 



(C4) 
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is the Dirac delta function extracting the transverse com- 
ponent. We obtain the equations of the radiation field as 



A = — n, 

1 



n 



where 



>cxT 



-X7 Z A 



P 



-(P + aA) 



•exTi 



(r) 



J dr 6 T (r-v')- J ex (r') 



(C5) 
(C6) 

(C7) 



is the transverse component of the current density (|A19|) . 
The equations of the polarizable oscillators are 



X = -(P + aA), 
P 



(C8) 



P = -p£j a X [aX] L / duj v^Q^j + aV^ ex , 

£ o P Jo 

(C9) 

and those of the reservoir oscillators are obtained as 

Y w = i(CL+^X), (CIO) 
P 

Q w = -pu; 2 Y u . 



(Cll) 

From Eqs. (C5J), ijCG]) . and fC8]) . we obtain the 
Maxwell wave equation for the vector potential 



V 2 A — ^-A 



pa 



aX 



MoJcxT, 



(C12) 



which has the transverse component of the excitonic cur- 
rent density compared to the same kind of equation in 
Ref. 0. Using a relation between the longitudinal com- 
ponents of excitonic variables 



•cxL 



(r) 



cxL 



(r) = e o V0 ex (r) 



(C13) 



and that for the polarizable oscillators (|B4|) , we can 
rewrite Eq. (|C12[) as a wave equation for electric field 



V x V x E + -^E = /j aX - /i J c 



On the other hand, from Eqs. (|C8[) and (|C9[) . we obtain 
the differential equation of the polarizable oscillators: 



pX+pd> 2 X = aA [aX] L +aV</>, 



1 f 00 
, N - / d„- < t Q < 
P Jo 



(CIS) 



2. Laplace transform 

Next, we rewrite the equations of motion derived in the 
previous section to those for the forward Laplace trans- 
form of the variables 

/>oo 

fl(p) = / dte- pt Q(t), (C16a) 
Jo 



and for the backward Laplace transform 



dt e~ pt n(-t). 



(C16b) 



From these motion equations, we derive the ones for 
positive-frequency Fourier transform 

n+(oj) = \n(-iuj + 5) + n(iuj + 6)] . (Ci7) 

Z7T 



From the forward Laplace transform of Eqs. (|C10|) and 
(IC11I) for the reservoir oscillators, we obtain 



QUp) 



— j^X(p) 



p 2 + UJ 2 



(C18) 



and, by using Eq. (|C8j) at t = 0, that for the polarizable 
oscillators (IC15|) becomes 



a 
P 



pA(p)-i[aX(p)] L + V0 ox (p)| 



1 

P JO 



du « w Q tl ,(p)+pX(0) + -P(0). (C19) 



Substituting the former into the latter and using the for- 
ward Laplace-transform of the electric field 

E(p) = -pA(p) + ±- [aX(p)] L -V0 ox (p)+A(O), (C20) 

we obtain the same equation as Eq. (20) of Ref. 0: 

X(p) = - -x(p)E(p) + %(p)(aA(0) + ppX(O) 



a 



P(0) 



dw 



p 2 + dj 2 



w 2 Y w (0) - ^CUO) 



where 



(C14) x(p) 



P + W 2 / dw 

P Jo 



p 2 + CiJ 2 



(C21) 



(C22) 



is the background susceptibility. 

Next, we derive the forward Laplace transform of 
Maxwell wave equation (|C14[) for the electric field as 



V x V x E(p) + ^E(p) - pap 2 X{p) 



(0) + 

- /LtoaX(O) — fioapX.(0). 



-U(O) 



r E(0) 



(C23) 



Here using the Eqs. dE}, ([C5|, ([C8|, (|CT2|) . and 

(|C13I) . we obtain the relations between the variables at 
t = as 



X(0) = -A(0) + -P(0), 
P P 



(C24) 
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E(0) = c 2 V x V x A(0) + — A(0) + — P(0) J cx (0), 

(C25) 

E(0) = --n(O) + - [aX(0)] L - V^ cx (0). (C26) 

Using these relations, we can rewrite the Maxwell wave 
equation (|C23|) to 



V x V x E(p) 



f(p)E(p) 



(C27) 



where e(p) = 1 + x(p) is the background dielectric func- 
tion, and the operator 

J(p) = — V x V x A(0) - e PX{p)A(0) + 11(0) 

1 - ^xGOl X(0) - [oX(0)] L - -px(p)p(o) 



6 f°° V 

—PXip) / 

a J p 2 + 



w 2 Y w (0)- -Q w (0) 

(C28) 

is the same one shown in Eq. (27) of Ref.0- This operator 
depends only on the background system variables and 
momenta of Ref. 0- 

On the other hand, the backward Laplace transform of 
the Maxwell wave equation is obtained as 



V x V x E(p) 



tKpMp) 



c 2 

where the operator on the RHS 

J(p) = — V x V x A(0) - e o px(p)A(0) - 11(0) 



(C29) 



a 



1 - e ^jP 2 x(p) X(0) + [aX(0)] L - -PX(P)P(0) 



a J Q p 2 + to 2 



a 

*rY w (0) + -QUO) 
P 



(C30) 



is also independent from the variables associated with 
cxcitons. 



3. Fourier transform 

From the forward and backward Laplace transforms 
(IC27I) and (|C29[) of the Maxwell wave equation, we obtain 
that for the positive-frequency Fourier component of the 
electric field as 



V x V x E+(r,w) 



-e bg (r,w)E+(r,w) 



ili u Jo(r,w) + J+ (r,w) 



(C31) 



where £bg(r, = e(r, — \lu + 8) is the background dielec- 
tric function. The noise current density operator in our 
system is written as 

\uj 2 

J (r,w) = J(r,w) jlm[eb g (r, w)] 

7TCT 

x /" dr' G*(r,r»- [J e x(r',iw + 5) - e V<Mr',0)] , 

(C32) 

where G(r, r, u) is Green's tensor satisfying Eq. ©, and 
J(r, w) = ^- [J(r, -iw + 5) + J(r, iw + (5)] 

-^Im[e bg (r,o,)] J dr' G*(r, r', lu) ■ J(r', k, + 6) 

(C33) 

is the same operator shown in Eq. (45) of Ref. satisfying 
the commutation relations |(7J). By using the relation 



J ox (r,cj) = -KjP ex (r,u/) 
we can rewrite Eq. (|C31[) to ([2"T]). 



(C34) 



APPENDIX D: EVALUATION OF 
COMMUTATORS 

Although we derive the Maxwell wave equation con- 
sidering the excitons in App. [Cj there exists a problem 
in deriving the commutation relations of noise current 
density J (r,o;) and in describing the optical processes 
of excitons. These are the subjects of this appendiex. 
In Sec. ID 1[ we derive the self-consistent equations for 
the Laplace-transformed operators, and we evaluate the 
commutators of the noise current density in Sec. ID 21 In 
Sec.lD3l we consider nonradiative relaxation of excitons. 



1. Laplace transform of self-consistent equations 

From matter Hamiltonian (jTTJ) and interaction Hamil- 
tonian (|15p . neglecting the radiation contribution of the 
current density AT 0X (r)A 2 (r)/2 under weak excitation, we 
obtain the Heisenberg equation of excitons as 



fc^t) - j dr [X;(r) ■ A(r,i) - ^(r)0(r,t)] 

dr» [gUtyd^t) + 9 ;(n)dl(n,t)] , (di) 



where <p(r) = <fibg( r ) + <^ex( r ) is the complete Coulomb 
potential. In this section, we neglect nonradiative relax- 
ation and assume lu ~ lu^. Using the relations between 
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X M (r), p p (r), and P M (r) (Eqs. flM6]) and ifXlS]) ) and the 
Laplace transform of the electric field 

E(r,p) = -pA(r,p)- V0(r,p) + A(r, 0), (D2a) 
E(r, p) = pA(r, p) - V0(r, p) - A(r, 0) , (D2b) 



the forward and backward Laplace transforms of Eq. (|D1|) 
are respectively derived as 

(hio^ — fko — iS) iu> + S) = —ihb^(0) 

+ J dr V*(t) ■ [E(r, -iw + S) - A(r, 0)] , (D3a) 

(hiu^ — frio + i5) 6 M (iw + S) = i^&p(0) 

+ J dr P* (r) • [E(r, ho + 5) + A(r, 0)] . (D3b) 

Adding these two equations, we obtain Eq. (|2"8f . the 
Fourier transform of the excitons' motion equation. On 
the other hand, using the relations 



eoV0ex(r) = PexL(r), 



(D4) 



and 



J ex (r,p) = pP cx (r,p) - P cx (r, 0), (D5a) 
J ox (r,p) = -pP ex (r,p) +P cx (r,0), (D5b) 



the Laplace transforms (|C27j) and (|C29j) of the Maxwell 
wave equation are rewritten as 



V x V x E(r, -ho + S) ^-ebgfr, w)E(r, -ho + 5) 

c z 

= ifi uj [J(r, -ho + 5) - P cx tO, 0)] 
+ /i w 2 P GX (r, -ko + S), 



(D6a) 



or 



V x V x E(r, ho + 6)- ^-e£ g (r, w)E(r, k> + 5) 

= i/i w [ J (r , ho + 5) + P oxT (r, 0)] 

+ /i w 2 P ox (r,iw + <5). (D6b) 

Here, under the RWA, the Laplace transforms of the po- 
larization density operator can be written as 

P ox (r, -ho + 5) ~ V^b^-ho + <*), (D7a) 
P ox (r, ho 

+ 6)~J2'P l *(r)h(^ + S ) ( D7b ) 



in the same manner as its Fourier transform ([29|) . Substi- 
tuting Maxwell wave equations (|D6|) in motion equations 
(|D3[1 of excitons, we obtain the self-consistent equation 
set for the Laplace transformed operators: 

^ [(ton - ^)<W' + •Aft,n'(^)] bf+'i-iuj + S) 

= ^(-^0 + 6), (D8a) 



= ^(iw + *), (D8b) 
where the operators on the RHS are 



P^iuo + S) = ifioco J dr I dr' V*{t) ■ G(r, r', w) 
• [J(r',-i^ + <5)-P cxT (r',0)] 
-i%(0)- f dr p;(r)-A(r,0), 



(D9a) 



^(iw + *)=iA*ow |dr / dr' 7>*(r) • G*(r, r', w) 
• [J(r',i^ + 5) + P cxT (r',0)] 
+ i%(0) + / dr P^(r) • A(r, 0). (D9b) 

In order to derive (|D8bj) . we used the representation 

A;,Ju) = -fiOJ 2 J dr J dr' P;(r).G*(r,r»-7V(r'), 

.( D1 °) 

which comes from definition (|35p and the reciprocity re- 
lation 



Gy(r,r» = G^r', r, w) 



(Dll) 



shown in Eq. (1.53) of Ref. 13. Adding self-consistent 
equation sets (|D8[) . we obtain the one for the Fourier 
transform of the exciton operators as 



^ [(^V - fko)S^/ + A^<(co)] = /3 M (u;), 



where the operator on the RHS is 



/^M = ttM- [uj + S ) + ttP^ + 5 ) 

ZlT Z7T 



(D12) 



+ ^ E [A^M - a;,Jco)] v(i W + 5). (Dis) 



Now, we verify that this operator is equivalent to the 
RHS of Eq. (1371). Using Eqs. (lC32ll . (lD4l) . (jD5l) . and 
(jTTJ) , the background electric field (124)) is rewritten as 



E+(r,w) = i/i w y dr' G(r, r', w) • J(r', w) 

-£/dr' [G(r,r»-G*(r,r»] 

• [w 2 P ex (r', iw + 5) + iwP cxT (r', 0)] . (D14) 
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Since, from Eqs. (IC33|) and (fTTj) , we can rewrite the first background field (ID14|) can be rewritten as 
term as 



J dr' G(r,r»-J(r» 

= dr' G(r,r , ,w).J(r / ,-iu> + 5) 

+ i- |dr'G*(r,r',c) ■ J(r', iu; + £), (D15) 



E+(r,w) 



2tt 



and the RHS of Eq. (J3TJ) becomes 
dr P*(r).E+(r,w) 



2tt 

2 
1 

+ 2^ 



dr' G(r, r', w) • [J(r', -iu + S) - P cxT (r', 0)] 
dr' G*(r,r»- [j(r', ice + 5) + P cxT (r', 0)] 
dr' [G(r, r', «) - G*(r, r', «)] • P cx (r', iw + 5), 

dr / dr' V;(r) ■ G(r, r', w ) • [J(r', -iu + 5) - P cxT (r', 0)] 



(D16) 



^ /dr/ dr ' ^(r).G'(r,r', W ). [J(r' ; ic + <5) + P cxT (r', 0)] 



(D17) 



r 



We can find that this is equivalent to (|D13[1 by expanding 
+ 5) and ^(iw + 5) with Eqs. (|D9)) . 



the first two terms and the following terms of (|D9b[) are 
commutable as 



2. Commutation relations 



In the representation of background field (|E)14[) or 
(|E)16|) . the information of excitons during t < is re- 
flected via operator P ox (r,kj + 5). In order to evalu- 
ate the commutation relations of noise current density 
Jo(r, cj), we first discuss the motion equations of exci- 
tons. 

The exciton motion during t < is described by 
Eq. (|D8b[) with (|D9b[) . Using W(w), the inverse matrix 
of coefficient matrix S(u) defined in Eq. (136|) . we rewrite 
Eq. (|D8b|) to 



J(r,iw + 5)+P cxT (r,0), 
i%(0) + y dr' p;(r')-A(r',0) 



0. (D21) 



Then, since P OX T(r, 0) and A(r, 0) are commutable with 
themselves, we obtain the expression to be evaluated: 

%{iu + S),0 v (iu' + *')}*] = h 2 S^ 

+ fi 2 ^' J dr J ds J ds ' / dr ' p p( r ) -G*(r,s,w) 

• [J(s, iu + 6),{J(s', iu' + 8')}*] ■ G(s', r', u') ■ T v {v'). 

(D22) 



b,{iu + W P» h> + *)• (° 18 ) Us i ng the relation shown in Eq. (Bl) of Ref. 



First of all, we evaluate the cummutation relations of 
$n(kj + 5). Since, from Eqs. (|C30p and (US]), we obtain 
the relations 

[J(r, iw + S), A(r', 0)] = iM T (r - r'), (D19) 



[%iw + 5),{J(s',iu/ + <5')} t ] 
= v / VV - IV 2 6(s - s' 

' 0h V~7Tx [ £ b g ( s ^) ~ £ b g (s,w')] lJ(s-s'), 



a; — a;' — i<5 



[P exT (r,0),6 /i (0)] = -[^(r)] t 



(D20) 



(D23) 
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we obtain 



where 



[iq'loJ J ds J ds' G*(r, s,w) 

• [J(s, ic^ + 6), { j(s', vJ + S')Y] • G(s', r', w') 



- u/ - i<5 



c^G*(r,r»-w' 2 G(r,r',w') . (D24) 



From this relation, Eq. (ID22|) is evaluated as 



[/3 M (ic + ( 5),{^K + <5')} t ] 



i<5 



(D26) 



We also obtain the relation 



because this commutator has a nonresonant denominator 
com pared to ([D26]) . Therefore, from Eqs. ([D18]) . (|D26|l . 
and (|D27[) , we obtain the commutation relations of the 
backward Laplace-transformed exciton operators 



[^(iw + <5),{V(iu;' + 5')}t] 

+ 5), V(iw' + 5')] = 0. (D28b) 



On the other hand, using relations (|D4[) and (|D5p . 
noise current density Jo(r,w) defined in Eq. (|C32j) is 
rewritten as 



J (r,w) = J(r, uj) - Im[eb g (r, u)] I dr' G*(r,r',wj 

7TCT ' 



P ex (r'> + 5) + -P oxT (r',0) 



(D29) 



Then, the commutator with its Hcrmitc conjugate is eval- 
uated as 



Xi = - Im[e bg (r,w)] 

x J ds G*(r,s,w)- [p^s.iw + ^.^rV)}*] , 

(D31) 



/3 



Im[e bg (r'V)] 



x / ds' 



J(r,^),{P ox (s'>' + 5)}t -G(s',r',c'), 

(D32) 



w 3 a/ 3 

X 3 = "^T Im [ £ bg(r,w)] Im[e bg (r',w' 
x /" ds J ds' G*(r,s,u/) 



P ex (s,icj + (5) + -P ex T(s,0), 



{P ex ( S ',iu/ + <5)}t--P exT (s',0) 



G(s',r>'). 

(D33) 



First, from Eqs. (|D9j) and (|D18j) . the commutator ap- 
pearing in Xi becomes 



[p ex (s,iu; + <5),{j(r',u/)}i 

= ^p M (s) w;,,» | ds' t>;,(s') 



i/xow / ds" G*(s',s",w) • 3(s",[lu + 5) + A(s',0), 



(D34) 



{j(r>')}t 



Since J(r',o/) is defined as Eq. (|C33j) . what we must 
evaluate is 



J(s",iu; + ( 5),{J(r',w')} t 

= ±- \J(s", iu + S), {J(r', -iu;' + 5) + J(r', iu;' + 5)}t] 

27T 

+ ^Im[e bg (r',u/)] | dr" 

[J(s", k, + 5), { J(r", iu/ + 5)}t] • G(r", r', u/). 

(D35) 



From Eqs. (Bl) and (B3) of Ref. 7, the first term is eval- 
uated as 



Jo(r,u,),{j (r',a/)}tj = [j(r, w ), {J(r' f w')}t 

+ X!+X 2 +X 3 , (D30) 



[J(s", vu + 8), {J(r', -io/ + <5) + J(r', iu/ + <5)} f ] 
2feoI m [ e bg(r', 



i(u> — lo' — iS) 



-lS(s"-r'). 



(D36) 
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Since the second term obeys Eq. (|D24|1 . we obtain 

ds" G*(s',s",w) • [j(s",\w + S),{3{r',w')}^ 
e film[e bg (r',w')]a/ 3 



-G(s',r>' N 



m lu(lo — w' — iS) 



(D37) 



From Eqs. (52) and (53a) of Ref. [j], the last term of 
Eq. (|D34|) becomes 



ds'P* (s')- A(s',0),{J(r',a/)}t 



hw 



J Im[ £bg (r>')] |ds' [p;,(8')] T -G(8',r', W '). 

(D38) 



Therefore, Xi is evaluated as 



Im[e bg (r,a;)] Im[e bg (r', u/) 



x^|ds G*(r,8, W ).7» P (8) 



x / ds' 



- w' - iS 



G(s',r',w'). 

(D39) 



From its complex conjugate, we obtain 



hw 2 w' 2 



Im[e bg (r,w)] Im[e bg (r', w')] 



— ■ — -\r OJ [T^{S T 

ui — uj + id 



x J ds G*(r,s,w) 
xVH / ^' P;,(s')-G(s',r', W ') 



(D40) 



Since, from Eqs. (fD7b|) . (|D9b)) . (fD18]) . and (fD28|) . the 
commutator appearing in X3 is written as 



P ex (s,kj + <5) + -P exT (s,0), 



{P ex (s',L/ + <5)}t-^P exT (s',0) 



fj,,fj,' 



w — w' — 10 a/ L F J 



w — w' — 10 a; 



we can find that 



Xi + x 2 + x 3 = 0. 



(D41) 



(D42) 



Therefore, we obtain commutator (|22a[) of the noise 
current density. From the same kind of calculation 
and neglecting the nonresonance terms, we also obtain 
Eq. ph) . 



3. With nonradiative relaxation 

In this section, we perform the same kind of calcula- 
tion as in the previous section considering a nonradiative 
relaxation process of excitons. Considering the reservoir 
oscillators, the Laplace transforms of Heisenberg equa- 
tion (ID1D of excitons are derived as 



(hw^ — hw — id) b/j,(— iw + 5) 

= -i%(0) + J dr V;{r) • [E(r, -iw + S) - A(r, 0)] 

- / dO [ 9/J (n)(J M (n, -iu + S)+ g*(tt)dUtt, -iu + 5)] 
Jo 

(D43a) 

(hujfj, — hw + iS) b^{iuj + 5) 

= i%(0) + /" dr V;{r) • [E(r, iw + 6) + A(r, 0)] 

p CO 

■ - / do [ 5/1 (n)d M (n, io; + 5) + (0)4(0, uj + ^)] . 

Jo 

(D43b) 

On the other hand, we obtain the motion equation of 
reservoir oscillators 

iftj^n,*) = nn d„(si,t) + g*{n) [& M (t) + 6+(*)] , 

(D44) 

and its Laplace transforms 
(/if2 — fiw — i<5)d M (0, — iw + 8) 

= -illd^n, 0) - g*{Sl) [b^-iw + 5) + 5J, (-iw + 5)] , 

(D45a) 

(Hfl -hw + i6)dn(Sl,iu + 8) 

= ihd^n, 0) - g ;(si) [b^ + s) + &t(i^ + <*)] . 

(D45b) 



Substituting Eqs. (|D45[) into (|D43[) and neglecting the 
nonresonant terms, we obtain 

[hw^ - hw - i7 M (w)/2] b^—iw + 8) 
= -i%(0) + J dr pj(r) • [E(r, -iu, + 5) - A(r, 0)] 
+ D fl (-iu + 6), (D46a) 

[Sw^ - tiw + i7*(w)/2] b^iLu + 8) 

= i%(0) + /" dr 7>;(r)- [E(r,iw + <5) + A(r,0)] 

+ £> M (iw + <S), (D46b) 
where the operators on the RHS are defined as 



D^-iw + tf) 



dO 



du(fi,0) 



l, O — uj — iS 

jg;(") 



dO „ 4( n »°)' ( D47a ) 
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D^iu + S) = - 



+ 



dQ 



Q — uj + iS 



i Q + u> — 16 
and the relaxation width is 

2 



774(0,0), (D47b) 



17m M _ 







|3m(")| : 



ftfi — — i5 Ml + htjj + iS 

(D48) 

Substituting Laplace transformed Maxwell wave equa- 
tions (|D6|) into Eqs. (|D46[) . the self-consistent equation 
sets for the Laplace transformed exciton operators are 
obtained as 



£S^ s ,(u;)V(-iu + £) 



= ^(-iu + 5) + + S)= (3^ s (-iuj + 6), 

(D49a) 

^£{S^(uj)}%,(ibJ + 8) 
m' 

= ^ (iw + <S) + £»^(i^ + 5) = /3f s (iw + 5) , (D49b) 
and the one for the Fourier transform is 

E ^M'M VM = + A» = /% b »> (D50) 

m' 

where operator V^uj) is defined as 

V^(uj) = i- + 5) + D^{\lo + 5)] 



1 i7 M H+i 7 ;(a;) 
2tt 2 



& M (iw + <5). (D51) 



On the other hand, Fourier transformed motion equation 
(|47| of the excitons is obtained by adding Eqs. (|D46p . 
Further, by substituting Eq. (|23|) into it, we can directly 
obtain self-consistent equations (|D50|) . 

From commutation relations (fll))) of cL(f2), the ones of 
the operators D^(— iui + 5) and D^iw + S) are calculated 
as 



[D lt (-m; + S),{D lt ,{-k/ + S)}i] 
= feW i7 M (a;) +i7*(u; / ) 
w - a/ + i<5 2 



cj — uj' — i<5 2 



[^(-iw + ^A*'^' + <*)}*] 

= ^M,M' 17m (^) ~ 17m (^0 

(w + itf) - (uj' + iS) 2 



(D52a) 



(D52b) 



(D52c) 



and neglecting the nonresonant terms, we obtain 
[D„(-iLj + 5),D^(-iu;' + 8)] 

= [D li (-h) + 6), + 5)] =0. (D52d) 

From these relations, we obtian the commutation rela- 
tions of $* hs (iuj + S): 



uj — uj' — 16 



(D53) 



uj — uj' — 16 2 
—J— - [SJ* (a/) - {^»}*] , (D54) 



and 



[ft hs (iuj + 5)J$ s (iuji + 5)]=Q. 



(D55) 



They have a good correspondece with Eqs. (|D26|) and 
(|D27[) . which are discussed without the nonradiative 
relaxation. Therefore, the commutation relations of 
b^iiuj + S) have the same form as that of Eq. (|D28[) ex- 
cept for replacing W^.^(uj) with Wjj~*,(ui), and relations 

([22]) of Jo(r,Ct>) are not changed even by considering the 
nonradiative relaxation. 

Next, we evaluate the commutators of X> M (w) defined 
in Eq. (|D51[) . Here, we apply the Markov approximation 
to the reservoir oscillators interacting with excitons, and 
assume that correction term 7^(w) is a real value. From 
Eqs. (|D52p and Dirac's equivalence (x ± it)) -1 = Pa; -1 =p 
mS(x), we find 



[D^iuj + 5) + D^iuj + S), 

{D^(-iuj' + 6)+D fl ,(iuj' + 6)Y] 
= ^m,m'^( w - uj'^irh-j^uj), 

[D^-iuj + 6) + D^{iuj + S), 
D M , {-iu' + S) + (iuj' + 5)] = 0. 

On the other hand, Eq. (|D49bj) gives the relation 

[D^-iu + S)+ D^iuj + 5), {b^{iuj' + 5)}t] 

7m ( w ) - 7m K) 



(D56a) 



(D56b) 



(uj + iS) ~ (uj' + i<5) 



ih 



7:(w)+7m(w') 



w; b ;,( w '). (D57) 



Then, we obtain Eq. (|48a| and also Eq. (|48bj) by neglect- 
ing the nonresonent terms. 

Finally, we verify the independence of source operators 
Jo(r,w) and T>^(uj). Since V^(uj) is defined as (|D51[) , we 
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obtain the expression of the commutator 
2^H,{J (iV)}t" 

D^-ko + 6) + D^iw + 5), {J (r, w')}t 



'M 

1 

2^ 



i 



7mM +7uH 



2tt 



6„(iw + <5),{J (r,a/)}+ 



(D58) 

Here, from Eqs. (|D29j) and (|D57j) . the first term becomes 



47T 2 C 



■ Im[e bg (r,w')] 

7m H - 7m( w ') , 7mH + 7m( w ') 



(w + id) - (w' + i8) to - lu' - iS 
£ W^"') / dr ' K'( r ') ■ G(iW). (D59) 



The commutator appearing in the second term of (|D58[> 
is written as 

b^iou + d),^,^] = [^(ia> + *),{J(r,<y)}+ 



Im[e bg (r,u/)] / dr' 



6 M (iu; + 5), 



{P cx (r', vJ + <5)}t - -^P cxT (r', 0) -G(r', r, w'). 

(D60) 

From Eq. (|D39)) , the first term on the RHS of (|D60| 

fiu/ 2 



Im^a/)]^^^)}* 



x / dr' 



— u>' — iS 



G(r',r,w'), 

(D61) 



and from Eqs. (|D9b|) , (|D49bp . and (|D28|) . the second 
term is calculated as 



hw' 2 



Im[e bg (r,u/)] 



m' 



+ {W^(cj)V [K'^t] ■ G ( r >^')- (D62) 
Therefore, commutator (|D60[) is evaluated as 

hu>' 3 Im[e bg (r, w')] 



- oj' - id 



[6 M (iu; + <5),{Jo(r,a/)} 1 



From this and Eq. (|D59j) . we obtain Eqs. (|49|) by as- 
suming 7^(w) = 7u(^) and neglecting the nonresonant 
terms. 



APPENDIX E: SELF-STANDING MODES AND 
RETARDED CORRELATION FUNCTION 



As discussed in sections [V] and I VII commutation re- 
lations of exciton operators in the Fourier transformed 
Heisenberg representation are written as Eq. (|42[) in 
terms of inversed matrix W abs (cj) of coefficient S abs (w) 
of the self-consistent equation set. As the analogue of 
the relation between the electric field and Green's ten- 
sor G(r, r', oj) satisfying Eq. ([5]), we can understand that 
the elements of W abs (lu) identify with the Fourier trans- 
forms of retarded correlation functions D^*, of exciton 
operators: 



Rex / 



dtD^,(t)e^, (El) 



wehre 







t > f 
t < f 



(E2) 



This identity indicates that W abs (w) has no pole in the 
upper half side of the complex w-plane. On the other 
hand, the poles {lu\} in the lower half w-plane charac- 
terize the self-standing modes of the exciton-polaritons 
satisfying 



det[S ab > A )] = 0, 



(E3) 



because it makes inverse matrix W aDS (w) singular at lu = 



APPENDIX F: EQUAL-TIME COMMUTATION 
RELATIONS 

Commutation relations of equal-time Heisenberg op- 
erators should keep the form of those of the Shrodinger 
operators. This means that the relations 



K{t),bUt) 



= s, : 



[P ex (r,i),P cx (r',t)] =0, 



[E(r,t),E(r',i)]=0 



(Fl) 



(F2) 



(F3) 



should be derived from the commutation relations of the 
Fourier transformed Heisenberg operators. Moreover, 
from Eqs. (|2"U|) and (|C5[) . the electric field is represented 
as 

e E(r,i) = -n(r,i)-e V<Kr,f), (F4) 

and since II(r,t) and A(r, t) satisfy commutation rela- 
tion (|C1|) . the relation 

[e E(r, t), A(r',t)} = iM T (r - r') (F5) 
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should also be derived. For local dielectric media, the 
same kind of calculation has been performed by Knoll, 
Scheel, and Welsch (KSW)^ 

From the time-representation of exciton operator 



Kit) 



duj b^{oS) e 



(F6) 



and relations (|52")) for w-representation, the equal-time 
commutation relation is written as 



h 

\2n 



Aw[W^)-{W^)Y 



(F7) 

In the limit \lo\ — ► oo, as indicated in App. A.l of KSW 
work^ it is known that £b g (r, ui) — > 1 and 



lim 



-G(r,r', W ) = -5(r-r'). 



(F8) 



Then, due to the orthogonality of ^(r) shown in 
Eq. ((27]) and the relation with LT splitting A£ T = 
|"PM| 2 / e bg(^v)eo, the limit of correction term ([35]) is 

lim A^(lu) = - fdr V;(r) ■ V^(r) (F9) 

|w|->oo £q J 



e bg (w /J )A^ T . 



(F10) 



On the other hand, the limit of nonradiative width 7 M (w) 
defined in Eq. (|D48[) is 7^(1^) — > 0. Therefore, coefficient 
matrix S abs (co>) becomes diagonal as 

lim S* h °, (w) = [hw^ + e bg K) A^ T -hu-iS\ 5 MtM , , 



|oj| — >oo 

and of course its inverse matrix is also diagonal: 



M 



lim <£M 



/kj^, + eb g (w Al )A^ T — ?uj — i<5 



(Fll) 
(F12) 



As mentioned in App. [El since W abs (w) has no pole m 
the upper half w-plane, the integration over the real axis 
is evaluated as 



a* ran 



— A 

h 



(F13) 



We can find that this equation reproduces commutation 
relation (JfTJ from Eq. (|F7|) . 

Next, we verify commutation relation (|F2[) of exci- 
tonic polarization P ex (r, i). Although we approximate 
its positive-frequency Fourier component as Eq. (|29[) in 
Sec. El here we describe it without the RWA as 

P+ (r, W ) = £ [p^r) 6» +7>;(r) {^(-w)}* 

(F14) 

This representation keeps the following relation derived 
from the definition of Fourier transform Eq. ([5]): 

P+(r, w) = P" (r, -co) = {P+(r, -«)}t (F15) 



From this relation, the time-representation can be writ- 
ten as 



P«*(r,f) 



(111.' 



P+ x (r, W ) e" 



P cx (r,w) e 



'not 



(F16) 



Pi(r,a;) e" 



dwP-(r,a;) e iwt . 

(F17) 



Therefore, the equal-time commutator becomes 

[P cx (r, i),P cx (r',i)] 

«U' / do/ e~ iut [p+ (r, w), P"(r', w' 



E[nwn( r ')-nwn(r')] 



(Fis) 



We cannot obtain the second term from w-representation 
(1291) with the RWA. Equation (|F18|) is also obtained using 
Eqs. (HSJ) and jFTJ) directly. From Eq. (|F18|) . we can 
reproduce equal-time commutation relation (jF2[) : 



[P cx (r,t),P cx (r',i)]^, 

= E[^( r ){^'( r ')r-c. 

A" 

= ^[<0|Pi(r)| M )(/i|P|;(rO|0) 

A" 

= (0| [i*(r),p£(r')] |0> = 0, 



(F19) 
(F20) 
(F21) 



where = x,y,z, \fi) = &TJ0), and |0) indicates the 
ground state of the medium. 

Next, we verify relation (|F3[) of the electric field oper- 
ator. From the representation of the positive-frequency 
Fourier component of electric field (|2"3")1 and that of exci- 
tonic polarization (jF14p . instead of Eq. (|39p. we can also 
write the electric field operator without the RWA: 

E+(r,u;)=E+(r,u/) 

+ 53[£ M (r, W )S p (w)-^(r,-a;){6 M (-a;)}t] J (F22) 



where we use the relation derived from Eqs. (14T))) and (|4"l"l) 
with G(r, r', w) = {G(r, r', -w*)}*: 



^(r,w)=^;(r,-w) 
^(r,w)=g;(r,-a;). 



(F23) 
(F24) 



Since the time-representation of the electric field is also 
written in the same form as Eq. (|F17|) . we can evaluate 
its equal-time commutator as 
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[E(r,i),E(r',i)] = / dec 
do; 



Ef(r,w),Eo(r / ,w' 



I 



(F25) 



The first term is the equal-time commutator between the 
electric fields in the background medium; then, it be- 
comes zero as indicated by SW or as calculated by KSW. 
Based on the fact that W^,(u>) and G(r,r',w) have no 
pole in the upper half w-plane as discussed in App. [Ej 
we evaluate the second term using the residue theorem. 
Since G(r, r',u>) becomes (|F8I) in the limit of \u>\ — > oo, 
from Eqs. (|40|) and (|4Tj) . we can obtain 



lim £ M (r,Ci») 

\uj\ — >oo 

lim ^(r,w) 



-r M (r)/eo, 



(F26) 
(F27) 



From these and Eq. (|F12j) . the second term of Eq. (|F25|) 

becomes 

^E / dcJ w) *0 M *V (r', w) 

^E^W^M- (F28) 

Calculating the other terms in the same way, we obtain 
[E(r,t),E(r',t)] 

= ^E tnW^;(r') - V;(r)r,(r')] = 0, (F29) 



then, equal-time commutation relation (|F3[) is repro- 
duced. 



Last, we verify relation (IF5|) . As discussed in App. A. 2 
of KSW work, we can write the time-representation of the 
vector potential as 



A(r,i) 








poo 




= lim 




1 ds 












— ioji 


w 


do; — 






-oo 


iw 




{■ oc 




= -p 


/ da; 





-E+(s,w) +H.c. 



5 T (s - r), 
(F30) 

dsE+(s,w)-<5 T (s-r), (F31) 



where P indicates the principal value integration. From 
this representation, the commutator between the electric 
field and the vector potential becomes 



[e E(r,i), A(r',i)] = -e„ / dwP/ dw' ( . . } / ds 



E+(r,^),E -(s,^') .^T(s-r') 



2tt 



dw 



ds 



S;(v,-uj) W^,(-lo) ^(s-^-F^r-w) W^, B^a,-^ ■* r (s-r'). 

I 



(F33) 



The first term is the same kind of commutator for back- 
ground medium, then it becomes ifi.<$x(r — r') as verified 
by KSW. Therefore, all we have to do is verify that the 
other terms become zero. First, in the limit of \u>\ — > oo, 
they becomes zero because of factor u -1 comparing to 



the last four terms of Eq. (|F25j) . In the limit of |u;| — > 0, 
due to the equation shown in KSW work, 



lim 



+o c 1 



J ds G(r,s,w) • «5 T (s-r') = 



(F34) 
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we can find that the terms also become zero: Therefore, Eq. (|F33|) reproduces equal-time commutation 

relation (lF5l). 



lim — / ds S T (s - r') • £Jb, lj) 

= lim / ds5 T (s-r')-^ r u(s,o;) = 0. (F35) 
\t»\->o c z J 
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